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We report the first fully self-consistent calculations of the nonlinear optical properties of super- 
lattices. The materials investigated are mono-layer superlattices with GaP grown on the the top 
of InP, A1P and GaAs (110) substrates. We use the full-potential linearized augmented plane wave 
method within the generalized gradient approximation to obtain the frequency dependent dielectric 
tensor and the second-harmonic-generation susceptibility. The effect of lattice relaxations on the 
linear optical properties are studied. Our calculations show that the major anisotropy in the optical 
properties is the result of strain in GaP. This anisotropy is maximum for the superlattice with max- 
imum lattice mismatch between the constituent materials. In order to differentiate the superlattice 
features from the bulk-like transitions an improvement over the existing effective medium model is 
proposed. The superlattice features are found to be more pronounced for the second-order than the 
linear optical response indicating the need for full supercell calculations in determining the correct 
second-order response. 



o 



S 

i 

-a 



in 



o 
o 



I. INTRODUCTION 



Semi-conducting strained superlattices (SLs) are potential materials for applications in optical communications 
involving switching, amplification and signal processing. In particular III-V semiconductor hetero-structures and SLs 
have attracted a great deal of interest mainly due to the possibility of tailoring band gaps and band structure s] 1 ! 2 ! 3 ! 4 ! 5 
by variation of simple parameters like superlattice period, growth direction and substrate material. With the de- 
velopment of new techniques like the strain induced lateral ordering process^ and existing methods like molecular 
beam epitaxy^ and low-pressure chemical- vapor deposition^ it is possible to grow and tailor these SLs. Thus a 
great deal of experimental work has been devoted to these materials. The unusual optical behaviour of AlP/GaP 
O \ SLs has been extensively studied using photoluminescence, magneto-photoluminescence)2ii2iii optical absorptionji* 
X-ray diffraction 1 ^ and refractive index measurements. 14 InP/GaP, being one of the material combinations which 
\ spontaneously constructs the SL under specific conditions, has been subject to numerous experimental works 1 ^* 1 ^ 1 ^ 
^ • including cathodoluminescence experimentsi&ASi 2 ^ and studies on the influence of pressure, SL period and barrier 
thicknesa 21 ! 22 ! 23 on the optical transitions. Likewise GaAs/GaP^ 4 * 2 ^ and GaAs/AlAo 26 i 27 i 28 i 29 i 30 i 31 i 32 have also been 
extensively investigated in the past. 

Much of the theoretical work done to explain these interesting physical properties of SLs has been largely concerned 
with the understanding of the electronic band structure. For example, the effect of strain on the band gap, the 
■ band offset problem and the possibilities of engineering it as well as the interface energy and band structure have 
been studiedi 3 ^* 3 ^ 3 ^* 3 ^ 3 ^* 3 ^ 3 ^ 4 ^ 4 ^ 4 ^ The linear optical properties of SLs have also been determined theoretically. 
Franceschetti et al. have calculated the pressure coefficients of optical transitions in InP/GaP SLs^ 3 - The pseudopo- 
tential method within the local density approximation (LDA) has been used by Kobayashi et al. iA to calculate the 
optical transition strengths of the (AlP)„/(GaP) m type multi-quantum-well SLs and for GaAs/AlAs by Yee et idJ± 
Confinement effects on optical transitions for various superlattice periods in GaAs/AlAs has been discussed by Schmid 
et aZ. 4 - The sp 3 s* tight-binding method has been used for the study of optical properties and the indirect-to-direct 
band gap transition of AlP/GaP^ 4 ^ and GaAs/AlAs SLsi 4 ^ Botti et nlM- have employed the linear combination of 
bulk bands method to determine the optical properties of GaAs/AlAs SLs and electroabsorption properties of these 
SLs have been studied by Kawashima et alM> Shibata et al£Z used the pseudopotential method to determine the 
. £h oscillator strength and band structures of the AlP/GaP SLs. 

Due to the simultaneous breaking of the time reversal and inversion symmetries at surfaces and interfaces, the 
$_i ' nonlinear optical properties are more sensitive than the linear optical properties and so the second harmonic generation 
(SHG) by some of these SLs has also been calculated and determined experimentally^ The major theoretical work 
in this direction was done by Ghahramani et a l£2£££££& They used the effective-medium-model (EMM) to determine 
the linear and nonlinear optical properties of (GaAs)„/(GaP) n , (Si) n /(Ge)„ and (GaAs)„/(AlAs) m SLs. They also 
employed the non-self-consistent linear-combination-of-Gaussian-orbitals (LCGO) method within the LDA to calculate 
the band structures and optical properties of these compounds. In these works Ghahramani et al. conclude that away 
from the absorption edge, most of the features in the linear as well as nonlinear optical properties are due to the bulk 
transitions which are well modeled by the EMM. 

There have been very few studies of the linear optical properties of SLs and even fewer calculations of the nonlinear 
properties. Also the existing theoretical work is neither self-consistent nor takes into account the effect of lattice 
relaxations on the SL optical properties. Hence, there is a need for a fully self-consistent calculation of the nonlinear 
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and linear optical properties of superlattices in conjunction with a study of the lattice relaxation effects. The aim 
of the present work is to calculate the linear and the nonlinear optical properties of III-V mono-layer SLs using the 
state-of-the-art full-potential linearized augmented plane wave method (FPLAPW). The materials investigated are 
(InP)i/(GaP)i, (AlP)i/(GaP)i and (GaAs)i/(GaP)i. The EMM is also tested for all these materials. For InP/GaP 
and AlP/GaP, this model fails to reproduce even the features away from the absorption edge. So the EMM is 
reanalyzed and a modification is proposed. 

The paper is arranged in the following manner. In Section II we present the details of the calculations. Sections 
IIIA and IIIB deal with the linear and second-order optical response of mono-layer SLs, respectively. Section IV 
provides a summary of our work. Appendix A lists all the formulae used for the calculation of the SHG susceptibility. 
Details of the modifications to the EMM are presented in Appendix B. 

II. METHODOLOGY 

Total energy calculations are performed using the FPLAPW method implemented in the WIEN2k codeS Thereby 
the scalar relativistic Kohn-Sham equations are solved in a self-consistent scheme. For the exchange-correlation 
potential we use the generalized gradient approximation (GGA) derived by Perdew and WangiS 

The detailed formalism for the determination of the linear dielectric tensor e(u>) = + ie2(w) within the 

FPLAPW formalism has been presented before^ The susceptibility for the second harmonic generation \^ (2w, u, uj) 
has been calculated using the an extension to this program^ The formulae for calculating x( 2 )(2w,w,w) have also 
been presented beforei 61 i 62 i 63 i 64 They are rewritten to improve the computational efficiency, and we compare the 
computing time for both sets of formulae in Appendix A. The SHG susceptibility obtained satisfies all the theoretical 
sum rules presented by Scandolo and BassaniS 

All the calculations are converged in terms of basis functions as well as in the size of the fc-point mesh representing 
the Brillouin zone. The linear optical properties are calculated on a mesh of 500 fc-points in the irreducible Brillouin 
zone (IBZ) and for the second-order susceptibility a mesh of 1500 fc-points in the IBZ is used. 

TABLE I: Superlattice lattice parameters in A units. a^=a x =a y and a±—a z which is the direction of crystal growth. 



Superlattice 




a± 


(InP)i/(GaP)i 


4.150 


5.660 


(AlP)i/(GaP)i 


3.933 


5.457 


(GaAs)i/(GaP)i 


3.998 


5.553 



In all the superlattices GaP is grown on top of other group III phosphates (InP and A1P) or GaAs as the substrate 
material. The superstructures are constructed in the (110) direction, which is equivalent to (100) direction in these 
materials. All these structures are treated in tetragonal unit cells. The details of this unit cell for such SLs has been 
discussed before (Fig. 1(a) in Ref. l6t|) . The lattice constants of the SL are estimated using the macroscopic elasticity 
£ neor j i 34 J 42 4 67 anc j are p resen ted in Table I. The experimental lattice constants of the constituent materials (InP, GaP, 
A1P and GaAs) are taken from Ref. |68|. 

III. RESULTS AND DISCUSSION 
A. LINEAR OPTICAL RESPONSE 

First, we shall discuss the importance of lattice relaxations and their effect on the linear optical properties. Total 
energy and force calculations performed for the unrelaxed atomic positions of the mono-layer InP/GaP and AlP/GaP 
SLs show considerable forces on the atoms, so the relaxation of the structure is needed. On complete relaxation there 
is a gain in energy of 0.195 and 1.124 eV per formula-unit for (InP)i/(GaP)i and (AlP)i/(GaP)i, respectively. For 
the case of (GaAs)i/(GaP)i the atomic relaxations were not needed. 

The imaginary part of the dielectric function, £2(w), for all the SLs under investigation are presented in Figure 1. 
All the optical results presented in this paper are scissors corrected^ using the scissors operator calculated by taking 
the difference between the theoretical and experimental band gap of the SL. Among the SLs under investigation for 
the relaxed structures, the anisotropy in e2(io) is maximum for (InP)i/(GaP)i (Fig. 1(a)), where the lattice mismatch 
between InP and GaP is 7.38 % , while for (AlP)i/(GaP)i (Fig. 1(b)) the anisotropy is minimum and so is the lattice 
mismatch (0.2%). (GaAs)i/(GaP)i (Fig. 1(c)) lies in between the two with a lattice mismatch between GaAs and 
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FIG. 1: Components of the imaginary part of the dielectric tensor for the relaxed (SLR) and un-relaxed (SLUR) superlattice 
for (a) (InP)i/(GaP)i, (b) (AlP)i/(GaP)i and (c) (GaAs)i/(GaP)i. The important features in ef*(a>) for the relaxed SLs are 
labeled with the uppercase letters. 



GaP of 3.66 %. These trends indicate that most of the anisotropy is due to the strain in the growth material rather 
than to the lower symmetry of the superlattice. These results are in accordance with similar observations previously 
made for GaAs/AlAs and Si/Ge SLsi££*£& On the other hand the anisotropy in ea(w) of the unrelaxed SLs shows 
a maximum for the SL formed from the materials with minimum lattice mismatch (AlP/GaP). Also £2(0;) of the 
unrelaxed SLs is lower in terms of peak heights compared to the relaxed SLs. Hence, the lattice relaxations are very 
important in determining the correct anisotropy and peak heights and thus all the further calculations are performed 
for the relaxed SLs. Our results for €2(0;) of (GaAs)i/(GaP)i are qualitatively in reasonable agreement with previous 
calculations by Ghahramani et al. 53. However, the magnitude of the response in the present results is nearly two times 
that obtained previously. This is expected since in their work Ghahramani et al. pointed out that the LCGO method 
is not able to reproduce the correct magnitude of the optical response and their results could be underestimated by 
a factor of 2. The other differences in the two theoretical results are as follows: (1) The present calculations give one 
sharp peak around 4.75 eV (labeled E), while previous results show a broad multiple peak structure in e| z (w) and 
ef a: (w) around the same energy. (2) In the energy region 3.2 - 4.6 eV ef^w) is greater in magnitude than e| z (w) in 
the present work, where as the previous results show the opposite. These differences can be due to the difference in 
the method of calculation. Ghahramani et al. have used a non-self-consistent LCGO method within the LDA, while, 
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the present calculations are fully self-consistent using the FPLAPW method within the GGA. 




FIG. 2: Frequency dependent obtained by the SL calculations compared to the EMM results for (a) (GaAs)i/(GaP)i, 

(b) (InP)i/(GaP)i and (c) (AlP)i/(GaP)i. In addition the data for the bulk GaAs (a), InP(b), AlP(c) and the strained bulk 
GaP ((a), (b) and (c)) are also presented. 

A macroscopic model for the dielectric function of SLs has been suggested by Ghahramani et alJ& In this model 
the SL is considered to be constructed of slabs of an unstrained bulk substrate material (InP, GaAs or A1P in the 
present case) and a strained bulk material (GaP in the present case) grown on top of the substrate. This model is 
called the effective medium model (EMM). We have used this model to estimate the macroscopic averaged response 
and the results are presented in the Figure 2. Like Ghahramani et alM we find that the EMM generates reasonably 
good intensities but poor peak positions in (GaAs) \j (GaP) i (Fig. 2(a)). In the optical spectrum of (InP)i/(GaP)i 
(Fig. 2(b)) and (AlP)i/(GaP)i (Fig. 2(c)) the peak positions are better but the intensities are less well described. 
Another point that should be noted is that for (AlP)i/(GaP)i, where the lattice relaxations are important, the EMM 
calculated response is closer to that of the unrelaxed SL than to the relaxed SL response, but the same is not true for 
(InP)i/(GaP)i such that no consistent picture emerges. 

In order to find the reason for this behaviour, €2 zz (uj) for bulk InP, GaAs and A1P and strained GaP are also 
presented in Figure 2. We find that if the major peak positions in the dielectric function of the constituent materials 
are well separated in energy the EMM gives two separate peaks. As can be seen for the (InP)i/(GaP)i the major 
peak in the optical spectra of InP and strained GaP are around 5.0 and 6.65 eV respectively. The EMM also gives 
two peaks of almost equal strength around 5.0 and 6.65 eV. The separation between the major optical peaks of GaAs 
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and strained GaP is 0.76 eV and between A1P and strained GaP it is 0.80 eV. The EMM in these cases results in 
broadened and thus weakened peaks. It is surprising that the SL response is seemingly dominated by the substrate 
material with the growth material having little ((GaAs)i/(GaP)i) to nearly no contribution ((InP)i/(GaP)i) (Fig. 
2). 




Energy hco (eV) 



FIG. 3: ef (oj) calculated using the full SL and the SCEMM and ef^w) calculated using the SCEMM for (a) (GaAs)i/(GaP)i, 
(b) (InP)i/(GaP)i and (c) (AlP)i/(GaP)i. In addition the data for the bulk GaAs (a), InP(b), AlP(c) and the strained bulk 
GaP ((a), (b) and (c)) are also presented. The important features in e|*(o>) for the relaxed SLs obtained by the SCEMM are 
labeled with the lowercase letters. 

The purpose of the model is to find the bulk-like features in order to facilitate the detection of true SL features 
by comparing the model results with the full SL calculations. But EMM takes into account the effect of junction 
formation only partially by using the strained lattice parameters for one of the constituent materials. This is useful 
for reproducing the anisotropy. The effect of the junction formation on the band gap, however, is ignored which is 
taken to be the experimental gap of the unstrained bulk. As the number of SL layers increases this effect diminishes 
but in the case of mono-layer SLs it is pronounced. So the failure of the model in certain cases is not a surprise. 
Nevertheless, in order to provide a simple model for predicting the bulk-like features in the SL optical properties 
on the basis of its constituent materials, the shortcomings of the EMM can be fixed in the following way: The 
calculations are performed first for the unstrained growth material (GaP in the present case) and the scissors operator 
is determined by comparison with the experimental data. Now assuming that GGA (or LDA) underestimates the 
gap consistently in the unstrained and strained bulk material this same scissors operator is used to correct the 
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calculated linear (and nonlinear) optical response of the strained growth material (GaP in the present case). The 
average of this response with the response of the bulk substrate material (A1P, InP or GaAs in the present case) 
is then taken to determine the effective optical response of the SL. We call this new model the strain-corrected 
effective medium model (SCEMM). The details of the SCEMM are presented in the Appendix B. The results of the 
SCEMM along with the dielectric function of the substrate and strained growth material, corrected by the scissors 
operator calculated as stated above, are presented in the Figure 3. As can be seen the model reproduces most 
of the features in the SL spectra of all the compounds under investigation. The remaining discrepancy in peak 
positions is due to the fact that the experimental SL gap is lower than the average gap of the constituent materials. 
More importantly the SCEMM shows that the SL response is not just dominated by the substrate material, as 
indicated by the EMM, but has features from both constituent materials. This fact is best seen in the case of 
(AlP)i/(GaP)i where the SCEMM peak labeled a in Fig. 3(c) is dominated by the strained bulk GaP transitions 
and c is bulk InP like. The anisotropy in 62(0;) determined by the SCEMM confirms our earlier observation that 
most of the anisotropy is due to strain in the growth material rather than to the lower symmetry of the superlattice. 
The features not reproduced by the model are referred to as the SL features in the present work. The following 
features can be identified SL effects: SX1(~ 2.5eV), SL2(~ 3.75eV), SX3(~ 4.25eV), SX4(~ 5.6eV), SX5(~ 6eV) and 
5L6(~ 7eV) in (InP)i/(GaP) 1 , 5L1(~ 2.85eV), SX2(~ 5.25eV), SL3(~ 5.8eV) and SL4(~ 6.5eV) in (AlP)i/(GaP)i 
and SX1(~ 1.88eV), SX2(~ 3.25eV), SX3(~ 5.25eV) and SL4(~ 6.57eV) in (GaAs) 1/ (GaP) 1 SLs in Fig. 1 

TABLE II: Static dielectric constants ef x (0) and ef z (0) and the components of the SHG susceptibility in the static limit 
Xl vx (0) and xf *(0) in xl(T 8 csu for the relaxed SLs compared to the SCEMM and the EMM results. 



Superlattice 


ef x (0) 
SL SCEMM EMM 


ef*(0) 
SL SCEMM EMM 


xT(o) 

SL SCEMM EMM 


x7 z (o) 

SL SCEMM EMM 


InP/GaP 


9.27 9.60 8.51 


9.00 9.30 8.17 


13.89 15.57 9.47 


10.42 15.55 9.29 


AlP/GaP 


8.65 8.25 8.27 


8.86 8.28 8.30 


6.34 6.08 6.14 


7.51 6.10 6.16 


GaAs/GaP 


10.48 9.76 9.23 


10.19 9.54 8.93 


16.65 11.4 9.14 


15.56 11.40 8.95 



The real part of the dielectric function in the static limit is a directly measurable quantity. Our calculated results 
of e xx (0) and £i z (0) for the relaxed SLs are presented in Table II. Two things should be noted. First the anisotropy in 
ei(0) follows the trend of the lattice mismatch with a maximum of 2.96% for (InP) 1 /(GaP) 1 and a minimum of 2.60% 
for (AlP)i/(GaP)i while (GaAs) 1/ (GaP) 1 lies in between the two with an anisotropy in ei(0) of 2.80%. Second, in 
case of (AlP)i/(GaP)i and (GaAs)i/(GaP)i the lattice constant for the substrate material (GaAs and A1P) is larger 
than that of the growth material (GaP) and so ef x (0) is greater than ef z (0), while it is just the opposite for AlP/GaP 
SL. The results of ei(0) using the SCEMM are closer to the full SL calculations (worst by a maximum of 6.8% ) than 
the results obtained using the EMM (worst by a maximum of 12.3%). A more important test for any model, however, 
comes in the determination of the nonlinear optical properties as they are much more sensitive to the small changes 
in the bandstructure than the linear optical spectra. 

B. SECOND-ORDER OPTICAL RESPONSE 

The magnitude of the SHG susceptibility x {2^, for the mono-layer SLs along with the results of the SCEMM 

and the EMM are presented in Figure 4. The peaks in x^ (2w, u, to) can be identified to be coming from 2io and/or 

uj resonances of the peaks in the linear dielectric function. Therefore this identification of the peaks in the SHG 

susceptibility is done from the respective linear optical spectra. The identified peaks are marked in Figure 4 and the 

nomenclature adopted is M[xui) + N(yu>), which indicates that the peak comes from an xui resonance of the peak 

M with the yco resonance of peak TV in the linear optical spectra. For example, for (InP)i/(GaP)i (Fig. 4(a)) the 

(2) 

hump just below 1 eV, labeled A[2u>) in the Xxyz(2ui, u), ui) component for the SL comes from the 2u> resonance of 
the peak labeled A in the linear optical spectra (Fig. 1(a)). The peak at 1.15 eV is a SL peak, labeled SLl{2u>), 
coming from the 2ui resonance of the peak labeled SL1 in the e 2 (o>) plot (Fig. 1(a)). Similarly all the other features 
have been identified and marked in Figure 4(a). As expected the pure SL peaks SLl(2u), SL2(2u), SL3(2ui) and 
SLl(u>) + SL5(2lu) are absent in the SCEMM results. The features coming from resonances of the SL and bulk-like 
peaks such as SLA(2w) + SLl(u) + C(2tu) and C(u>) + SLl(2u) are underestimated by the SCEMM as features c(2u) 
and c(w), while the peaks coming from the 2lo and/or u> resonances of the bulk peaks are well reproduced by the 
SCEMM. The peak labeled a(2u) is more pronounced than the corresponding SL peak A(2u>) in accordance with a 
(Fig. 3(b)) being larger than A (Fig. 1(a)) in the linear optical spectra. In contrast, the peak labeled B(2u>) is smaller 
than the corresponding SCEMM peak b(2ui) + a(u) since the SCEMM peak is due to the ui resonance of peak a with 
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Energy hco (eV) 



FIG. 4: 4(a) The two inequivalent components of susceptibility for the SHG \xiyz(2io,uj,ijj)\ (thick solid line) and 
[Xi5L(2ai,u;, w)| (dashed line) obtained from the full SL calculations, along with the SCEMM (thick dash-dot line) and the 
EMM (dotted line) results of |x^( 2w . w » w )l for (InP)i/(GaP)i . 4(b) The same for (AlP)i/(GaP)i . 4(c) The same for 
(GaAs)i/(GaP)i. 



the 2u> resonance of the peak labeled b in the linear dielectric function (Fig. 3(b)). A similar feature assignment for 
the (AlP)i/(GaP)i and (GaAs)i/(GaP)i mono-layer SLs is also performed and the details are marked in Figs. 4(b) 
and 4(c). The same trends are observed in all the materials under investigation, with the SCEMM reproducing all 
the peaks other than pure SL features. This indicates that the SCEMM can be used to identify the structure in the 
optical spectra coming from the bulk-like transitions and hence facilitating the determination of the SL effects like 
symmetry lowering. The EMM on the other hand is not able to reproduce all the non-SL features and is therefore 
not as good a model for the mono-layer SLs. 

Our results for the magnitude of \^ (2w, u>, uj) of (GaAs)i/(GaP)i (Fig. 4(c)) differ from previous calculations^ 
(1) X^xyz oj, uj) is generally in good agreement with the previous work in terms of peak heights as well as positions. 
The only difference is the appearance of two extra features in the present work, i.e., a hump at ~1.25 eV (A(2oj)) and 
a peak at 1.75 eV (SL1(uj)). (2) In the present work the anisotropy in the nonlinear optical response is found to be 
more pronounced compared to the work of Ghahramani et al. As is well known, the nonlinear optical properties are 
more sensitive to the small changes in the bandstructure than the linear ones so any anisotropy in the linear optical 
response is expected to be enhanced in the nonlinear spectra. This is in accordance with our findings, while the 
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previous work shows the opposite. (3) Our results for x' 2 ' ) 2 y X (2w, w, w) differ substantially from the previous results, 
but since a detailed identification of the peaks is not presented previously it is difficult to point out the discrepancies. 
As mentioned earlier, these differences could be due to the different methods used. 

Finally the SHG susceptibility in the static limit is presented in Table II. The maximum anisotropy in X (0) is m 
(InP)i/(GaP)i (28.5%) and the minimum is found for (GaAs)i/(GaP)i (12.8%). (AlP)i/(GaP)i has an anisotropy 
of 16.8%. Neither EMM nor SCEMM give results close to the SL results. This is in contradiction with the previous 
findings, 53 " where Ghahramani et al. report the EMM to give results within 5% of the SL data. The failure of the 
SCEMM and the EMM in the present case is not surprising since both the models fail to reproduce some of the main 
features in the optical spectra leading to a substantial difference in the static limit. But still the SCEMM results are 
an improvement over the EMM. 



IV. SUMMARY 



To summarize, we have performed calculations for the linear and nonlinear optical properties of some III-V mono- 
layer SLs. We conclude the following: (I) The lattice relaxations play an important role in the determination of the 
correct anisotropy and peak heights of the linear optical spectra. (2) The effective medium model (EMM) does not 
reproduce the correct bulk features mainly because it assumes that the optical band gap of the material does not 
change with strain and is taken as the experimental direct gap of the unstrained material. (3) An improvement over 
this model is the strain-corrected effective medium model (SCEMM), which correctly averages the optical properties 
of the constituent materials of the SL and facilitates the identification of the SL features in the optical spectra. (4) 
The SCEMM results confirm that most of the anisotropy in e(u>), in the SLs under investigation, comes from the strain 
in the growth material rather than the symmetry lowering due to SL formation. (5) The anisotropy in 62 (w) and 
in particular in £2(0) follow the trend in the lattice mismatch being maximum for (InP)i/(GaP)i and minimum for 
(AlP)i/(GaP)i. (6) The small SL features in £2(0;) are greatly enhanced in the nonlinear spectra and the same is true 
for the anisotropy. (7) The SCEMM is able to correctly reproduce all the bulk-like features in the SHG susceptibility, 
whereas pure SL features are absent and features coming from the resonance of the SL and bulk-like transitions are 
underestimated. This leads to a failure of the SCEMM (as well as EMM) in the determination of X (0)- -But the 
SCEMM is still an improvement over the EMM. (8) The SCEMM is good for determining the linear optical properties 
of the SLs, but for the correct determination of the nonlinear optical properties SL calculations arc essential. The 
effects of strain due to junction formation and the lattice relaxations are expected to decrease with increase in the 
layer thickness. It would be interesting to compare two models in such a case and to find the SL period for which these 
effects become insignificant. (9) Although no experimental measurements of the SHG susceptibilty for the compounds 
under investigation exist, a detailed comparison of future experimental data with theoretical results would help in the 
identification of various features in the optical spectra, in particular, to highlight the effect of the interface formation 
in the SLs. This would lead to a better understanding of the physical properties which is one of the most essential 
ingredient for tailorable materials of technological importance. 



V. APPENDIX A: FORMALISM FOR THE SECOND-ORDER RESPONSE 



The formulae for the total susceptibility for the second harmonic generation (SHG) for clean semi-conductors have 
been presented before&S We note that Eq. B3 of Rcf 61 is incorrect and the correct form of this equation can be 
obtained from the sum of Eqs. B16b and B17 of Refl62t The susceptibility for SHG can be divided into three major 
contributions: the intcrband transitions Xinter(2w, u>, to), the intraband transitions Xintra(2w, uj, uj) and the modulation 
of interband terms by intraband terms Xmod(2u;, oj, uj): 



abc fn. , , , . ,\ e / r nn A r ml r ln\ f 2/„ m fml fin \ /-. s 

1,-^1 47T \Ulln — U) m l J l{U mn -2oj) (uj m i - Uj) [U) ln - OJ) J 



ri ni 



47T 3 



win i, * nm 

.nml 



L0 mn r nm{ r ml r tn 1 



fn 



fit 



Si 



fn 



,} 



^ fn m r^ m {r b ml rf n }(oj ml - u ln ) 



nm 



nml 



(2) 



9 



"~2 ? V { W ni r im{ r mn r 'nl} — ^ImX^X { r lm r mn } } 



.nml 



(3) 



for all rt ^ m ^ /. The symbols arc defined as 

A« m (k) = <„(k) - «« m (k) (4) 
with u° m being the a component of the electron velocity given as 

<mW=™» m (k)CW (5) 

and 

K TO (k)r^(k)} = ^« m (k)r^(k) + r* m (k)r° ,(k)) (6) 

The position matrix elements between states n and m, r" m (k), are calculated from the momentum matrix element 
p nm using the relation^ 

CJk) = . Km(k) ,, (7) 
zmw„ m (k) 

for all w„(k) ^ w m (k) and r° m (k) = otherwise, with the energy difference between the states n and m given by 

Tiuw = h(uj n - uj m ). (8) 

For the sake of clarity the k dependence of ui nm , p nm , V^ rn and A„ m in the Eqs. (l)-(3) are suppressed. Equations 
(1) to (3) are computationally very demanding since n, I and m run over all the bands. But this requirement can be 
relaxed by adjusting the equations slightly using the fact that the equations are symmetric in n, I and m and hence 
the indices can be interchanged. At the same time one can get rid of the Fermi functions. Eqs. (l)-(3) then read 



Xintcr 



nml k 



nm I ml In J 



i 



' Iml mn' nil 



ml)(Wn 
r nli r ln 



2w) 



i) (w; m — UJmn) 



(9) 



e 3 I ' I 



- 8 »E 2 ttt 9 / n m {^f n } + 2 2J -" mL 2 m<, /" J " T "'o T"^ ( 10 ) 



^ „a f„b J 1 i o r nm{ r ml r ln}( U 'rnl ^ln) 



e 1 

2?l " k l™!^^™"^ 



r a f b Ac \ 

Enm X mn mn J / 1 i \ 

7X77 77) < 



where f2 is the unit cell volume, Wk is the weight of k point and n denotes the valence states, m the conduction 
states and I denotes all states (I ^ m,n). We have used the Eqs. (9)-(ll) to calculate the total susceptibility. We 
also confirm that the real and imaginary parts of this susceptibility satisfy all the sum rules presented by Scandolo 
et alr& Some results of the benchmark are presented in the following table. All the calculations are performed with 
a mesh of 1000 points in the energy interval 0.0 - 4.0 eV on a SGI R12000 processor. 
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TABLE III: Number of fc-points in the IBZ, number of valence bands (VB), number of conduction bands (CB), CPU time (in 
minutes) needed to calculate the SHG susceptibility using Eqs. (l)-(3) and Eqs. (9)-(ll). 



k points 
IBZ 


Number of bands 
VB CB 


CPU time 
Eqs. (l)-(3) Eqs. (9)-(ll) 


385 


9 11 


11.02 2.67 


385 


9 5 


4.50 1.10 


146 


9 11 


3.50 0.88 


146 


24 33 


88.01 19.47 



IS 



VI. APPENDIX B: FORMALISM FOR THE STRAIN-CORRECTED EFFECTIVE MEDIUM MODEL 

The expression for calculating the xx component of the linear optical response using the effective medium model 

55 



4' V A) = i(el>, 5) + e|>, A)) 



(12) 



which can be obtained from the continuity of the tangential electric field across the interface Et = E^ and Et = E^ 
and the expression for the average induced polarization P e tf = \{P A + P B ). The longitudinal electric field on the 
other hand is given by c^jjEl = e A z Ef — e g z E B , and so the expression for e e JJ is given by 



24*(w,5)e| B (w,A) 



(w, 6) 



(13) 



Here it is assumed that A is the substrate material and B is the material grown on top of A and hence has the 
strained lattice parameters. e A x (ui,5) is the xx component of the frequency dependent dielectric tensor for material 
A (or strained material B) corrected with the scissors operator S (or A). Where 6 and A are calculated as: 



5 = XDG - TOGx 



(14) 



A = XDG - TOG s 



(15) 



where XDG is the experimental direct band gap and TOG is the optical bang gap calculated theoretically. The 
subscripts X and S indicate that the experimental lattice parameters and the strained lattice parameters, respectively, 
are used for calculating the theoretical optical band gap. The EMM assumes that even after straining the material B 
has the same gap as measured experimentally for the unstrained material. These expression are slightly altered for 
the SCEMM reading: 



eff 



(16) 



c eff 



M) 



2e A x (to,8 A )e B x (u J ,S B ) 
(u;,6 A ) + ef(uj,S B ) 



(17) 



where the symbols have the usual meaning. One extra set of calculations is needed for the unstrained B material to 
determine the optical band gap and hence the scissors operator. Then the calculations are performed for the strained 
material B and the band gap is corrected with the "old" scissors operator. This effective zz component of the dielectric 
tensor can be used to calculate the xyz component of the SHG susceptibility as follows: 
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